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Abstract 

Boson boson scattering and Higgs production in boson boson fusion will be actively investigated 
at the LHC. We have performed a parton level study of all processes of the type qiq2 — > q3q'iq5Q6l~^l~ 
using for the first time a full fledged six fermion Monte Carlo event generator which employs exact 
matrix elements at 0{a^^). We have examined Higgs production in vector boson fusion followed 
by the decay chain H — > ZZ including exactly all electroweak irreducible backgrounds. 

In the high mass region we have compared the case of a relatively light Higgs with the no-Higgs 
results. The integrated cross section for the latter case is more than twice that in the former for a 
minimum invariant mass of the ZV pair of about 800 GeV. We find, in a preliminary analysis at 
parton level that, summing up the muon and the electron channels, about 25 events are expected 
in the light Higgs case for L=100 fb^^. 
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I. INTRODUCTION 



The Standard Model (SM) provides the simplest and most economical explanation of 
Electro- Weak Symmetry Breaking (EWSB). Detailed reviews and extensive bibliographies 
can be found in Refs.jl, 2, ^, 4|. The only missing ingredient is the Higgs boson. The fit of 
the SM to precision EW data currently gives an upper limit on the Higgs mass of about 200 
GeV [5] while direct searches have established a 95% CL lower bound M(H)> 114.4 GeV Q]. 

In the SM the Higgs is essential to the renormalizability of the theory and is also crucial 
to ensure that perturbative unitarity bounds are not violated in high energy reactions. 
Scattering processes between longitudinally polarized vector bosons (V^) are particularly 
sensitive in this regard. Without a Higgs the V^^'s interact strongly at high energy, violating 
perturbative unitarity at about one TeV [Z]. If, on the contrary, a relatively light Higgs 
exists then they are weakly coupled at all energies. In the strong scattering case one is led 
to expect the presence of resonances in Vl^l interactions. Unfortunately the mass, spin 
and even number of these resonances are not uniquely determined . If a Higgs particle 
is discovered it will nonetheless be necessary to verify that indeed longitudinally polarized 
vector bosons are weakly coupled at high energy by studying boson boson scattering in full 
detail. 

At the LHC no beam of on shell EW bosons will be available. Incoming quarks will 
emit spacelike virtual bosons which will then scatter among themselves and finally decay. 
These processes have been scrutinized since a long time, going from the pioneering works 



m 



10 



papers m 



hich address boson boson scattering on a general ground, to the more recent 
0, 3 focused on the extraction of signals of vector boson scattering at the LHC. 
All previous studies of boson boson scattering at high energy hadron colliders, with the 
exception, to our knowledge, of |l^ and [l^, have resorted to some approximation, either 



the Equivalent Vector Boson Approximation (EVBA) [lu\, or a production times decay 
approach, supplementing a calculation of 

qiq2 q3q4ViV2 (1) 

processes with the, typically on shell, decay of the two vector bosons. There are however is- 
sues that cannot be tackled without a full six fermion calculation like exact spin correlations 
between the decays of different heavy particles, the effect of the non resonant background. 
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the relevance of the offshellness of boson decays, the question of interferences between dif- 
ferent subamphtudes. Without a complete calculation it will be impossible to determine the 
accuracy of approximate results. In Ref. [l4| this issue was discussed at lenghth, showing 
differences of the order of 10-20% in some important regions of phase space. The reliability 
of the EVBA approximation in the context of vector boson scattering has been critically 
examined in 



18|. It 



Recently PHASE a full fledged six fermion Monte Carlo has become available 
describes at 0(aem); using exact tree level matrix elements, all processes of the form PP 



qiq2 - 
LHC 
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QaQ^Qp.Ii' (where qi stands for a generic (anti)quark) which can take place at the 



19|]. The range of interesting reactions is however much wider. Processes in which 



both vector bosons decay leptonically have been extensively studied both for Higgs detection 
and for boson boson scattering and top physics. Besides, in order to obtain a full coverage of 
all semileptonic processes it is necessary to include all reactions with a charged lepton pair in 
the final state. This has required the calculation of additional amplitudes and an extensive 
improvement of the routines which pilot the integration and the generation of unweighted 
events. The result is a new code called PHANTOM 20|] which, at present, includes all processes 
with six fermions in the final state at 0(ag^) 



PP qiq2 /l/2/3/4/5/6- 



(2) 



The accuracy of tree level calculations can be sensibly improved. PHASE is being contin- 
ued also in this direction. In particular, electroweak corrections have proved to be sizable 



expecially for processes involving the Higgs boson, see for instance Ref. 



21| . A new code 



PHAST_NLO 22] will address 0{aem) electroweak radiative effects in six fermion physics. Both 



PHANTOM and PHAST_NLO are based on the methods of Refs. 



23 



24l | and adopt the iterative- 



adaptive multichannel strategy developed in 18| . In the following we present results obtained 
with PHANTOM. 

PHANTOM is an example of dedicated event generator which describes a predefined set of 
reactions striving for maximum speed and efficiency. Other recent examples of dedicated 
programs for LHC physics are Alpgen j25| and Toprex 26] . The complementary approach 
is given by multi-purpose programs for the automatic generation of any user-specified parton 
evel process. The following codes for multi-parton production are available: Amegic-Sherpa 
27|, CompHEP 28|, Grace-GrOppa 29], Madevent Q, Phegas & Helac 31], O'Mega & 



3 



Whizard |32|. 

Recently the 0(ag„as) QCD corrections to qiq2 — > q'3q'4e+z/e/i~z7^ a have been computed 
in [3^ and those to qiq2 q3qAe~^ e~ fi~^ fi~ and qiq2 —>■ qsq^e^e^i^fjMfj. in 34|- They turn out 
to be modest, changing the total cross sections by less than 10%. The smallness of QCD 
corrections is related to color conservation which forbids gluon exchange between quark lines, 
up to higly suppressed contributions generated by the exchange of identical quarks. The 
same apphes to QCD corrections to qiq2 — > q^qAQblel^l' which however do include a larger 
number of terms which could result in larger corrections. It should be pointed out that QCD 
effects are not limited to NLO corrections. In the context of vector boson pair production 
PP WW, gluon initiated 0{al) processes gg — * WW, proceeding through a quark loop, 
have been shown to be a sizable correction to the tree level reactions qq — ^ WW when 



realistic Higgs search cuts are imposed 



35| . On the other hand QCD WWjj production 



represents only a small background to Higgs searches in Vector Boson Fusion 36 1. 

Since in addition to VV scattering many other subprocesses are in general present in the 
full set of diagrams, as partially shown in Figs. [THSl it is not a trivial task to separate boson 
boson scattering from the EW irreducible background. In practice one has to deal also 
with other types of background to which QCD interactions contribute, but which however 
do not include any boson boson scattering term. We will refer to these processes as QCD 
background even though in general they will be a mixture of QCD and EW interactions. In 
this paper we are neglecting QCD backgrounds. It is clear that obtaining a good signal over 
EW irreducible background ratio is a prerequisite to any attempt at dealing with the QCD 
one. 

We are aware that much is still needed to obtain a thoroughly realistic assessment of the 
observability of these processes. Only a complete study including 4j/+Z~ at 0(«em'^s) 
0{al^a^) together with full detector simulation will be able to say the final word. In the 
meantime it is important that the tools available for simulation are sharpened as much as 
possible and that the viability of such analyses is demonstrated at 0(ag^) which includes 
all signal contributions. In our opinion, the results presented in the following and the event 
generator used to produce them represent a step forward in this direction. 
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FIG. 1: Vector boson fusion processes. 




FIG. 2: Examples of EW irreducible background to vector boson scattering processes. 

II. CLASSIFICATION AND CALCULATION 

For a complete analysis one needs to include all processes which contribute to final states 
with one charged lepton pair. Taking into account one lepton type, charge conjugation and 
the symmetry between the first and second quark families, the number of reactions can be 
reduced to 135. A given reaction, its charge-conjugate, the ones related by family exchange 
and those obtained with the product of the two transformations can be indeed described by 
the same matrix element; they differ by the convolution with Parton Distribution Functions 
(PDF). All processes which share the same total particle content, with all eight partons 
taken to be outgoing, can be described by a single master amplitude. As a consequence, all 




FIG. 3: Higgs boson production via vector boson fusion and Higgsstrahlung. 
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Group 


Type 


diagrams 


Group 


Type 


diagrams 


uu cc bb l^l^ 


4Z 


615 


ddssbbl^l^ 


4Z 


615 


uu uu bbl^l^ 


4Z 


1230 


uu uu uul^l^ 


4Z 


3474 


uu uussl~l~^ 


4Z 


1158 


uu uu ccl^l^ 


4Z 


1158 


uu bb bbl^l^ 


4Z 


1606 


uu cc ddl~l^ 


4Z+2W2Z 


821 


uu ss ssl~l~^ 


4Z 


1158 


uu uu ddl~l~^ 


4Z+2W2Z 


2126 


uu ssbbl~l~^ 


4Z 


615 


uu dd bbl~l~^ 


4Z+2W2Z 


880 


ddddddl-l+ 


4Z 


3474 


uu ddssl^l^ 


4Z+2W2Z 


821 


bbbbbbl-l+ 


4Z 


7506 


uu dd ddl^l^ 


4Z+2W2Z 


2126 


ddddssl^l^ 


4Z 


1158 


udscuul^l^ 


2W2Z 


484 


ddddbbl-l+ 


4Z 


1230 


ud scddl^l^ 


2W2Z 


484 


ddbbbbl-l+ 


4Z 


1606 


ud scbbl^l^ 


2W2Z 


265 



TABLE I: Classification of pp qq' ^ 4g + l^l^ processes. The first column shows the group 
list, the second the process type as discussed in the text, the third the corresponding number of 
diagrams. 

reactions can be classified into 22 groups which are enumerated in Tab. [H By selecting two 
initial quarks in each particle group, one obtains all possible processes. 

The calculation can be further simplified examining more closely the full set of Feynman 
diagrams. In some processes, fermions can be paired only into neutral currents (4Z), while 
in other cases they can form two charged and two neutral currents (2Z2W). Mixed processes 
are described by a combination of the two sets (2Z2W+4Z). 
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The 4Z amplitude was not previously available. The three basic topologies in which the 
Feynman diagrams appearing in the 4Z amplitude can be classified are shown in Fig. |H The 
numbers N4/N3/N2/N1/N0 under each topology indicate the number of Feynman diagrams 
described by this topology if 4/3/2/1/0 fermion pairs are taken to be massive, namely to 
have non zero interaction with the Higgs boson. The fiavour of all four fermion pairs are 
assumed to be different. The complete set of diagrams is eventually generated by exchange 
of identical particles. 

Each rectangle on both sides of the central boson topology is the sum of several subdia- 
grams representing all possible decays of an off shell 7 or Higgs boson to four outgoing 
particles as shown in Fig. O These sets of subdiagrams are evaluated only once, with a 
substantial efficiency gain, and then combined together in the end. 



III. PHYSICAL PROCESSES 



Boson boson scattering and Higgs production in boson boson fusion produce intermediate 
states with two bosons and two quarks as shown in Fig. [1] In this study we have only 
considered final states in which one Z boson decays leptonically to yU"^/^" and the other (either 
Z or W) hadronically. If both bosons decay hadronically the signal cannot be distinguished 
from the QCD background whose cross section is much larger. Final states where both 
vectors decay leptonically have a smaller rate and have been left for future studies. 

A number of event samples representative of all possible processes of the form PP — >■ 
have been produced with PHANTOM. In order to comply with typical 
acceptance and trigger requirements, the cuts in Tab. [TTlhave been applied. The acceptance 
cuts are standard. The wide separation in rj between the two tag quarks is a well established 
technique for separating the scattering signal from the background ll|, ll2|, ll3|. We have 
imposed a minimum invariant mass cut of 60 GeV on jet pairs rather than requiring a 
minimum AR separation. It is well known that at large pt the two jets from the hadronic 
decay of a W or Z boson tend to coalesce. This issue has been discussed, for instance, 
at length in the ATLAS TDR Isl] in connection with the observability of a heavy Higgs 
decaying to a WW pair, where a number of jet finding schemes have been studied. We defer 
to Sub-Sect. IIIIBI a discussion of the effects of a separation in Ai? for the class of processes 
under consideration. It should be kept in mind that selection procedures can be optimized 
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FIG. 4: Fundamental topologies associated with 4Z processes. Total number of diagrams 
2263/1251/803/615/579 
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FIG. 5: Decays of the off shell neutral bosons Z, 7, H 



in different ways for different analyses and that such optimization has to be performed at 
hadron level rather than at the more primitive parton stage we are discussing here. 

In most of this paper, since we are mainly concerned with extracting signals of boson 
boson scattering from the irreducible background, we select events using flavour information, 
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Pt (lepton)> 20 GeV 
|?7(lepton)| < 3 
E(quark)> 30 GeV 
PT (quark) > 20 GeV 
|r/(quark)| < 5 
M(/+/-)> 20 GeV 
M(qq)> 60 GeV 
A7?(tag-quarks)> 3.8 

TABLE II: Standard acceptance cuts applied in all results. Any pair of colored fermions must have 
mass larger than 60 GeV. Ar/(tag-quarks) is the absolute value of the difference in pseudorapidity 
between the two tag quarks. 



which will be unavailable in actual analyses, for the identification of vector bosons and top 
quarks. In Sub-Sect. IIIIBI we show that our results are not substantially modified if one 
adopts a more realistic procedure entirely based on invariant mass cuts. 
We have used the CTEQ5L [38] PDF set with scale 

Q'=M'^+lj:pi,. (3) 

where pTi denotes the transverse momentum of the i-th final state particle. 

Many subprocesses (i.e. ZW -^ZW, WW ^ZZ, ZZ ^ZZ, qb ^qtV) will in general 
contribute to a specific six fermion reaction, tt processes will not contribute to the 4ql^l~ 
channels but single top production with an additional neutral boson emission will be present. 

It is impossible to separate and compute individually the cross section due to a single 
subprocess, since there are large interference effects between the different contributions. We 
can however select all complete 2^6 processes which include a specific set of subdiagrams. 
For instance, ZW ^ZW with on shell bosons is described by four Feynman diagrams. These 
same diagrams, with the two incoming external vector bosons connected to the initial fermion 
lines and the two final ones connected to their decay products, constitute the ZW — >ZW 
set of 2^6 diagrams. Several sets can contribute to a single process and therefore the same 
process can appear in different groups; for example ud udfi^fi^bb will be included in 
both the WW — i>ZZ and ZZ ^ZZ groups. As a consequence of this multiple counting, the 

9 



Invariant mass of the two central quarks | 



10-' 



0) 1U 
Q. 



' 10"^ 



10' 




200 400 600 800 1000 1200 1400 

M(q,q',^*^) [GeV] 
Invariant mass of the two central quarks + |i" | 




200 400 600 800 1000 1200 1400 



M(q,q',^>) [GeV] 



FIG. 6: Invariant mass distribution of the two charged leptons and the two most central quarks, 
for different sets of processes. The upper plot includes the set of cuts described in the text. In the 
lower plot a further cut for vetoing top production is applied. 
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total cross section is smaller than the sum of the cross sections for the various groups. The 
upper part of Fig. [6] shows the invariant mass distribution of the two most central quarks 
(when ordered in pseudorapidity 77) and of the two leptons for all reactions which contain 
the different subprocesses as well as the distribution for the complete set of processes. We 
assumed M(H)=150 GeV. Notice that the Higgs peak is present in the ZW — *ZW curve. 
This is due to processes, like for instance ud udddfi^^~ which in addition to the ZW 
^ZW set of diagrams (e.g. ud — > ddWZ) include also diagrams describing Higgs production 
in s-channel (e.g. ud — > udZZ). 

The group comprising top diagrams has a large cross section. The lower part of Fig. [6] 
shows the same distributions after top subtraction. Top candidates are identified requiring 
a b-quark and two other quarks in the final state of the right flavour combination to be 
produced in a W decay, with a total invariant mass between 160 and 190 GeV. 

If no Higgs is present, all SM scattering processes between on shell weak vector bosons 
grow linearly with the center of mass energy squared, with the exception of ZZ ^ZZ which 
in this case is zero. This behaviour is in agreement with the low energy theorem (LET) 41|. 
The lower part of Fig. [6] shows that the ZZ ^ZZ component is relatively small compared 
with the total distribution and as a consequence does not represent too serious a background 
to searches for new physics signals. The group including the ZW — s>ZW set of diagrams gives 
the largest contribution. 

In Tab. IIIII we present the total cross section for pp —>■ qq' — 4g + fi~^mu~ with the 
standard acceptance cuts in Tab. [Ill for different Higgs masses. In Tab. IIVI we show the scale 
dependence of the total cross section for two Higgs scenarios, adopting the usual recipe of 
varying the scale by a factor of two in either direction. This leads to a 4% difference with 
respect to the central value. 



A. Higgs production and its complete EW background in PHANTOM 

PHANTOM is capable of simulating Higgs production in VV fusion together with all its EW 
irreducible background for all channels and for any Higgs mass and may be particularly 
useful in the intermediate mass range, below the WW threshold, where the production 
times decay approach cannot be used. Its improved treatment of the EW sector needs to be 
complemented by an accurate description of QCD dominated backgrounds and of the effects 
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Mh 


Cross section (fb) 


150 GeV 


7.33±0.02 


200 GeV 


13.67±0.03 


500 GeV 


9.89±0.02 


NoHiggs 


7.34±0.02 



TABLE III: Total cross section for pp qq' — > 4g + fi^mu with the standard acceptance cuts in 
Tab. HI] for different Higgs masses. 



Mh 


Q/2 


Q 


2Q 


200 GeV 


14.21 ±0.02 


13.67±0.03 


13.11 ±0.02 


NoHiggs 


7.57 ±0.02 


7.34±0.02 


7.01 ±0.02 



TABLE IV: Scale dependence of the total cross section for pp — > qq' 4q + fj.~^mu with the 
standard acceptance cuts in Tab. HH for two Higgs scenarios. Q is defined in Eq.(l3|) 



of QCD NLO corrections. 

Higgs production in VV fusion followed by Higgs decay to WW or ZZ is the second 
most abundant production channel over almost the full range of Higgs masses which will be 
explored at the LHC. It is regarded as the channel with the highest statistical significance 



for an intermediate mass Higgs 36 



391] . Previous analyses have focused mainly on the WW 



channel. For an intermediate mass Higgs the dilepton final state H — >■ I^I^(*) — > lulu is 
slightly favoured with respect to the H — > iyiy(*) — lujj channel because of the W + nj 
background which affects the latter. In the first case the main background comes from ti 
production followed in importance by EW WWjj production which is estimated to be about 
10% of the signal. In the second case the main background comes from W + nj followed 
by ti and EW WWjj production. QCD WWjj production can be reduced to be of the 
same order of magnitude as the EW contribution using acceptance cuts and can be rendered 
essentially negligible by a central jet veto which is not so effective in the EW case |l3l. |36 |. 

The production channel qq qqH, H ZZ l^l^jj has been examined in j40| while 
the channels qq qqH, H ZZ l^l^uu and ZZ l^l^l^l^ have been considered 
in Q. The m-jj and in^m- channels are particularly interesting because they allow a 
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direct reconstruction of the Higgs mass which in lujj final states must be extracted from 
the transverse mass distribution. 

As an illustration of PHANTOM capabilities, the four body invariant mass distribution of 
the /^"^/i^ pair and the two most central quarks in Aq^^fi^ final states in the neighborhood 
of the Higgs peak is shown in Fig. [7] for M(H) = 150 GeV and M(H)=200 GeV. The plot on 
the left is obtained by zooming into the area around the Higgs peak in Fig. [61 Both plots 
show the results for all reactions which contain the different subprocesses, as described in 
Sect. Illlt as well as the total distribution. Assuming a mass resolution of ±10 GeV around 
the peak, the EW irreducible background amounts to about 3%(13%) for M(H)=150(200) 
GeV. An order of magnitude assessment of the statistical significance of such a cross section 
for Higgs discovery can be obtained by comparison with the jjlu channel. The main reducible 
background is QCD V + nj production. Assuming that the effect of acceptance cuts are 
similar in the WW and in the ZZ channel, one can estimate the ratio of the significancies S 
in the two cases as 

S{l+l-jj) ^ cj,,H ■ BR{H ^ ZZ) ■ BRjZZ l+l-jj) ^ 
Sili^Jj) ^ cr^oH ■ BR{H ^ WW) ■ BR{WW lujj) ^oF^ ^ 

Since ai^^j/ai+i-^j ^10^ and BR{ZZ l+l-jj)/BR{WW lujj) ^ 1/3 we are left 

with 

S{in-jj) BR{H ^ ZZ) 



S{lujj) ' BR{H WW) 
For M(H)>200 GeV the ratio of the two branching ratios is about 0.5 and on the basis 



of the studies of the ATLAS 



36| and CMS 39|] collaborations for the WW channel, one 



expects a good statistical significance, of order five, for qqH^ H ZZ, ZZ l^l^jj- This 
naive estimate is in rough agreement with the analysis of 4^ which obtains significancies 
slightly below four in the mass range 200 < M{H) < 300 GeV. 

Below the ZZ threshold, the Higgs branching ratio to ZZ reaches about 0.08 at M{H) ^ 
150 GeV. For a luminosity of 30 fb^^ about 40 events are expected in the H i^l^jj, I = 
e, IX channel. Only a complete analysis including all backgrounds and full detector simulation 
could tell whether this is enough for establishing a Higgs signal in this range of masses in 
the vector fusion channel. 



13 



5 10-= 



10^ 



Invariant mass of the two central quarks + | 



■ Total 

■ WZ -> WZ 

zz -> zz 

■ WW -> zz 

Top 



110 120 130 140 150 160 170 180 190 
M(q,q',ti*n-) [GeV] 



10"' 



10' 



10-^ 



10- 



10' 



Invariant mass ot the two central quarks + | 



Total 

WZ -> WZ 

zz > zz 
WW -> zz 

Top 



160 180 200 220 

M{q,q',n*ti) [GeV] 



FIG. 7: Distribution of the invariant mass M{VZ) of the two candidate vector bosons for a Higgs 
mass of 150 GeV and 200 GeV. The contribution of the various subprocesses is evaluated as in 
Sect, imi 



B. The high VV mass region 

In the absence of firm predictions in the strong scattering regime, trying to gauge the 
possibilities of discovering signals of new physics at the LHC requires the somewhat arbitrary 
definition of a model of VlV/, scattering beyond the boundaries of the SM. Some of these 
models predict the formation of spectacular resonances which will be easily detected. For 
some other set of parameters in the models only rather small effects are expected, see for 
instance 

The simplest approach is to consider the SM in the presence of a very heavy Higgs. While 
this entails the violation of perturbative unitarity, the linear rise of the cross section with 
the invariant mass squared in the hard VV scattering will be masked by the decrease of the 
parton luminosities at large momentum fractions and, as a consequence, will be particularly 
challenging to detect. At the LHC, the offshellness of the incoming vector bosons will 
further increase the difference between the expectations based on the behaviour of on shell 
VV scattering and the actual results. For >10 TeV, all Born diagrams with Higgs 
propagators become completely negligible in the Unitary gauge, and the predictions for all 
processes in Eq.([2]) reduce to those in the Mh oo limit. In this section we compare this 
minimalistic definition of physics beyond the Standard Model, which we call the no Higgs 
scenario, with the predictions of the SM with a light Higgs. 
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FIG. 8: Invariant mass distribution for M(yZ) > 800 GeV. The full and long-dashed line refer to 
the no-Higgs case, the short-dashed and dot-dashed ones to M(H)=200 GeV. All results satisfy 
the standard acceptance cuts. The full and short-dashed curves present the results for our signal 
definition. For the long-dashed and dot-dashed histograms we have further required |ry(Z;;)| < 2 
and \ri{qv)\ < 2. 



An analysis of selection cuts capable to increase the difference between the no-Higgs and 
light Higgs cases could provide some guidance for the search of signals of new physics in 
boson boson scattering. 

As already mentioned, in the absence of the Higgs, all SM scattering processes between 
on shell weak vector bosons grow linearly with the center of mass energy squared, with the 
exception of ZZ —>-ZZ. Therefore all possible reactions, ZW — >ZW, WW — >ZZ, ZZ — »-ZZ, 
should be carefully investigated. 

An interesting possibility is to investigate whether there exist or not an elementary Higgs 
boson by measuring the VV cross section at large M(VV). Previous studies jl4| have shown 
that kinematical distributions are quite insensitive to the value of the Higgs mass provided 
it is much smaller than the invariant mass of the VV system. 

In order to isolate the VV fusion signal, and more generally two vector boson production, 
from all other six fermion final state processes and investigate EWSB, different kinematical 
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FIG. 9: Pseudorapidity distribution of the lepton pair for the no-Higgs case (fuh) and for M(H)=200 
GeV (dashed). All events satisfy M{ZV) > 800 GeV. 



cuts have been applied to the simulated events. 

First of all, single top production is vetoed as discussed in Sect. IIIII Second, the invariant 
mass of the two charged leptons has to reconstruct the mass of a Z, and is required to be 
in the range Mz ± 10 GeV. In VV fusion an additional W or Z decaying hadronically 
is expected to be present. Therefore events are required to contain two quarks with the 
correct flavour content to be produced in W or Z decay, with an invariant mass of ± 10 
GeV around the central value of the appropriate EW boson. If more than one combination 
of two quarks satisfies these requirements, the one closest to the corresponding central mass 
value is selected. This combination will in the following be assumed to originate from the 
decay of an EW vector boson. 

In a third step, in order to reject events which can be identified with the production of 
three vector bosons, the flavour content and the invariant mass of the two remaining quarks 
is compared with a W and a Z. If compatible within 10 GeV with either, the event is 
rejected. The events satisfying all these constraints will constitute the "signal" sample. 

In Fig.[8]we present the invariant mass distribution of the two charged leptons and the two 
jets associated with the vector boson decay for M(H)=200 GeV and for the no-Higgs case. 
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16 


(7,9) 


5 


(3,2) 


3.20 


1.2 TeV 


13 


(6,7) 


3 


(2,1) 


4.33 


1.3 TeV 


11 


(5,6) 


2 


(1,1) 


5.50 


1.4 TeV 


9 


(4,5) 


2 


(1,1) 


4.50 






Ai? = 


0.4 






Mcut 


NoHiggs 


M(H): 


=200 GeV 


Ratio 


800 GeV 


18 


(8,10) 


10 


(6,4) 


1.80 


900 GeV 


12 


(5,7) 


6 


(4,2) 


2.00 


1.0 TeV 


8 


(4,4) 


4 


(2,2) 


2.00 


1.1 TeV 


5 


(2,3) 


3 


(2,1) 


1.60 






Aii = 


0.5 






Mcut 


NoHiggs 


M(H): 


=200 GeV 


Ratio 


800 GeV 


12 


(5,7) 


8 


(5,3) 


1.50 


900 GeV 


8 


(4,4) 


5 


(3,2) 


1.60 


1.0 TeV 


5 


(2,3) 


3 


(2,1) 


1.60 



TABLE V: Number of events as a function of the minumum invariant mass of the ZV fJ-^f^^jj 
pair for L=100 fb^^. All events satisfy \r]{Zii)\ < 2 and |f?(q'v)| < 2. In brackets we show the 
contribution of the (ZW,ZZ) final states. 

A number of selection cuts have been studied in order to increase the difference between 
the two Higgs hypotheses. Simple requirements of centrality of the lepton pair and of the 
candidate second vector boson have proved to be the most effective. The pseudorapidity 
distribution of the charged lepton pair in the two cases is shown in Fig. [91 The long-dashed 
dot-dashed distributions in Fig. [8] have been obtained with the additional contraints that 
\ri{Zii)\ < 2 and |77((3'y)| < 2, where qy refers to the quarks which are associated with the 
vector boson decay. The corresponding distributions for the ZW and ZZ final states are 
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presented in Fig. [TOl The cross section for qqZW and qqZZ production are similar, however 
the discrepancy between the no-Higgs case and the M(H)=200 GeV is larger for the qqZZ 
final state. 
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FIG. 10: M{ZZ) (left) and M{WZ) (right) invariant mass distribution for M{VZ) > 800 GeV. 
The full and long-dashed lines refer to the no-Higgs case, the short-dashed and dot-dashed one 
to M(H)=200 GeV. The full and short-dashed curves present the results for our signal definition. 
For the long-dashed dot-dashed histograms we have further required \ri{Zii)\ < 2 and |r/((j'y)| < 2. 



In the upper part of Tab. |V]we present the number of events as a function of the minimum 
invariant mass of the n~^fJ'~jj system for L=100 fb~^ with the set of cuts shown in Tab. [Tll 
In brackets we also give the separate results for the ZW and ZZ final states. The number of 
events is smaller than the expected yield in the 4g/iz/ channel but the differences between 
the two Higgs hyphotheses are larger. In fact, similar ratios are obtained with comparable 
number of events. 

In Tab. |V]we also show the effect of requiring a minimum AR separation among colored 
partons. The expected number of events decreases sharply, by about a factor of two for 
AR = 0.4 and by about a factor of three for AR = 0.5, in the NoHiggs case. The statistics 
for a light Higgs is less affected because the vector boson distribution is less central in this 
case and it is precisely the jets originating from the most central and higher pt W's and 
Z's which are most likely to merge into one jet under the effect of a AR cut. As already 
mentioned in Sect. IIIIl this issue is well known and various approaches have been tried in the 
literature. v^TLAS [Sj (Sec. 9.3.1.3, 19.2.10.2), favors using a small cone AR = 0.2 for the 
determination of the jet centroid and then a larger cone AR = 0.4 for collecting the energy 
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flow of tfie jets. In QCD studies at tlie LHC a typical separation AR = 0.5 is adopted. 
In Ref. ^, [s^] a different approach has been proposed. First, jets with a total invariant 
mass in the neighborhood of the electroweak vector meson mass are selected. Then these 
jets are forced to divide into two subjets, by varying the separation parameter y in the k± 



scheme. For jets originating from a vector meson decay the scale yyp'^ at which the heavy 
jet splits into two subjets is typically close to the vector mass, while for standard QCD jets 
the splitting scale is much smaller. The subject of jet reconstruction algorithms is still lively 
debated. Since EW vector bosons are crucial to many investigations at the LHC, we expect 
that a scheme which allows to distinguish jets produced in the decay of high pt W's and Z's 
will be devised for this kind of specialized studies. 




0.06 

0.05 

„ 0.04 
a 

0.03 

D 
T3 

0.02 
0.01 



WZ final states 



noh 
mh200 



-2 

1z 



FIG. 11: Pseudorapidity distribution of the lepton pair in qqZZ (left) and qqZW (right) final 
states. The full line refers to the no-Higgs case, the dashed one to M(H)=200 GeV. In all cases 
M{VZ) > 800 GeV. 



At the LHC, the expected mass resolution for vector bosons decaying to dijets is about 
5 ^ 10 GeV, depending on the boson transverse momentum j^. This makes it quite difficult 
to separate ZZ f^~^fJ'~jj from ZW fJ'~^f^~jj final states on the basis of the invariant 
mass of the jet pair. It is therefore of interest to explore alternative means of separating the 
two final states. The low energy theorem 4l|] predicts that A{W~^W'' ZZ) = s/v"^ where 
A is the scattering amplitude, s is the usual Mandelstam variable and v is the coupling 
strength of the gauge current to the Goldstone boson, which in the SM coincides with the 
vacuum expectation value of the Higgs field. As a consequence of crossing symmetry then 
A{W Z WZ) = t/v'^. Therefore, averaging between the two orientation of the incoming 
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ZW system, and neglecting vector boson masses, one expects, in the reaction center of mass, 
da/dcos9 ~ const for W+W' ZZ and da/dcos6 ~ (1 + cos6^) for WZ WZ where 
6 is the scattering angle in the boson boson center of mass. This has led us to study the 
pseudorapidity distribution of the lepton pair in qqZZ and qqZW final states separately, as 
shown in Fig. [TTJ Despite the fact that we are not in the center of mass of the VZ system, 
that the incoming vector bosons are not on their mass shell and that we are plotting the 
distribution of an angular variable which is not the cosine of the center of mass scattering 
angle, the general prediction that the two final states have different distributions is verified. 
In the NoHiggs case, the ZW channel distribution is almost fiat in the region |?7| < 2 with 
small peaks in the forward and backward direction, as suggested by the LET, while the 
ZZ one peaks at zero rapidity. It is somewhat surprising, but quite welcome, that the two 
distributions are markedly different also when a light Higgs boson is present in the spectrum 
as expected in the SM. It should be mentioned that the angular distribution depends on the 
vector boson pair invariant mass. The plot for M{VZ) > 300 GeV, which we do not show, 
displays a similar, rather central, behaviour for the two processes. 

The selection procedure employed for Fig. [8] and Tab. |V]is not fully realistic: no flavour 
information will be available for light quarks and 6's will be tagged only in the central part of 
the detector. At this stage we want to isolate as much as possible the VV fusion signal from 
all other production channels, with the aim to identify a possible signal definition which could 
play the role which was played by CC03 at LEP2 [42], capturing the essence of the physical 
phenomena under investigation and allowing comparisons between the results from different 
collaborations. It is however of interest to investigate whether the relevant distributions 
are sensitive to the details of the selection cuts. In Fig. [12] we compare the invariant mass 
distribution of the two charged leptons and the two quarks associated with the vector boson 
obtained with two different methods. In the first case (dotted line) we select only quark 
pairs that have the correct flavour content to be produced in W or Z decay while in the 
second (dot-dashed line) we consider all quark pairs. In the more realistic setting the top 
veto is applied to any triplet of quarks with a total invariant mass between 160 and 190 GeV; 
moreover alia events in which two quark pairs with mass between M\y — 10 GeV and Mz + 10 
GeV are present are discarded. In both cases we identify the candidate vector boson with the 
quark pair whose mass is closest to the nominal vector boson mass. The two distributions 
differ by about 20% at small invariant masses but agree quite nicely at invariant masses 
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above 800 GeV, showing that our results based on Monte Carlo level flavour information are 
not seriously degraded when selection procedures closer to the actual experimental practice 
are adopted. For comparison we also present the invariant mass distribution obtained by 
identifying the two most central jets as the vector boson decay products before (solid line) 
and after (dotted line) top vetoing. 

~ VZ invariant mass | 

Totai(V=2q central) 




^" 200 400 600 800 1000 1200 1400 

[GeV] 

FIG. 12: Invariant mass distribution of the lepton pair and the two jets from boson decay for the 
no-Higgs case. The solid(dashed) line is obtained identifying the two most central jets as the vector 
boson decay products before(after) top vetoing. The dotted line is obtained requiring the correct 
flavour content for the jets identified as decay products of both the vector boson and the top. The 
dot-dashed lines is produced using solely invariant mass informations to identify the vector boson 
and the top decay products. 

IV. CONCLUSIONS 

In this paper we have studied all qiq2 — > QsqiQbQe^^ ^~ processes at order 0{a^^) at the 
LHC using for the first time a full fledged six fermion Monte Carlo event generator. We 
have studied Higgs production and its irreducible EW background in vector boson fusion 
followed by the decay chain H ZZ l^l^jj-, including exactly all electroweak irreducible 
backgrounds and in particular the interference of EW ZZ + 2j and ZW + 2j production with 
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the signal. The EW irreducible background in the neighborhood of the Higgs peak amounts 
to about 1.5%(6%) for M(H)=150(200) GeV. We have examined how simple kinematical cuts 
can be applied at generator level to extract the VV signal from the irreducible background. 
In the high mass region we have compared the case of a relatively light Higgs with the no- 
Higgs case. The integrated cross section for the latter is more than twice that in the former 
for a minimum invariant mass of the ZV pair of about 800 GeV. Summing up the muon 
and the electron channels, about 25 events are expected in the light Higgs case for L=100 
]h^^ . These results are encouraging and show that a more complete analysis, including all 
QCD backgrounds and an accurate study of jet separation algorithms in the high invariant 
mass region, is worthwhile. 
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